Glass-like two-level systems in minimally disordered mixed crystals 
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THz spectroscopy is used to identify a broad distribution of two-level systems, characteristic of 
glasses, in the substitutional monatomic mixed crystal systems, Bai_2:Caa;F2 and Pbi_a;Cai,F2. In 
these minimally disordered systems, two-level behavior begins at a specific CaF2 concentration. The 
concentration dependence, successfully modeled using the statistics of the impurity distribution in 
the lattice, points to a collective dopant tunneling mechanism. 

PACS numbers: 63.50.-l-x, 78.30.Ly, 61.43.Fs 
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The low temperature thermal properties of glasses and 
certain mixed crystals are dominated by a low lying 
distribution of excitations usually described within the 
phenomenological tunneling model as two-level systems 
(TLS) P|. Although this TLS model is used to interpret 
the physics of disordered systems; justifications for either 
its assumptions or the well-known quantitative univer- 
sality of fundamental physical characteristics have not 
been found 0, A related soft-potential model may 
have a greater range of applicability 0, \^ , but remains 
completely phenomenological. Other microscopic mod- 
els which emphasize the potential energy landscape 
have not yet provided new insights into the low tempera- 
ture anomalies. After the discovery of TLS in chemically 
disordered mixed crystals 111 it became apparent that 
progressively disordered crystals, in which, unlike glasses, 
the saturation of TLS behavior had not been reached, of- 
fered good opportunities for exploring the fundamental 
mechanisms responsible for the broad distribution of low 
energy states So far TLS behavior in crystals has 
been demonstrated only in systems for which disorder is 
produced either by a low symmetry component with ori- 
entational degrees of freedom, or by chemical disorder. 
Extensive work on fluorite mixed crystals shows broad 
TLS-like behavior for aliovalent dopants in the range of 1 
- 45 mol % P, 13 . There the substitutional dopant cation 
introduces charge compensating interstitial anions while 
at the same time producing local strain and electric fields 
due to the different components. 

In this Letter we report on a THz study of the low 
temperature properties of two isovalent mixed crystal 
series where glass-like TLS have been observed. Both 
Tia,i-xCa,xF2 and Pbi_a;Caa;F2 have the fluorite structure 



and with some difficulty, single crystals of high quality 
can still be produced for x < 0.1. The appearance of TLS 
is consistent with the ideas that (1) the local dopant con- 
centration varies statistically and (2) a critical concentra- 
tion is required. 

Due to the decomposition of both Bai_:rCa:rF2 and 
Pbi_a;Ca2,F2 solid solutions under cooling, a high cool- 
ing rate is used during crystal growth. Solid solutions 
Bai_:rCa:rF2 (x = 0.02, 0.04, 0.06, 0.08), crystal diam- 
eter 22 mm, crystal length 30 mm and Pbi_a;Ca2:F2 {x 
= 0.02, 0.04, 0.06, 0.08, 0.10), crystal diameter 10 mm, 
crystal length 30 mm were grown by the Stockbarger- 
Bridgeman technique. The melt was maintained for 1.5 
hours at T = 1400 °C (Bai_^Ca:,F2) and 2 hours at T 
— 950 °C (Pbi_i,Ca2,F2). The crystals were grown at a 
rate of 10 mm per hour at T = 1400 °C (Bai-^^Caj^Fz) 
and T = 950 °C (Pbi_j;Ca:rF2). After the crucible 
was transported over the gradient zone the furnace was 
switched off and the crystals were cooled to room tem- 
perature during 8 hours for Bai_a;Caa;F2 and 4 hours 
for Pbi_2;Caa;F2. These rates are faster than the typical 
crystal growing cooling times of 16-24 hours. The result- 
ing samples, which have large single crystal domains, are 
heavily strained, but optically transparent. 

The dopant concentration was measured by powder x- 
ray diffraction of end chips as well as by Raman scat- 
tering along the entire length of each sample. X-ray 
diffraction measures the average lattice constant, and 
Raman scattering measures the shift of the T2g symme- 
try phonon mode. Figure ^ shows the Raman and x-ray 
data converted into the mol % CaF^ concentration, lin- 
early extrapolated via Vegard's law [lj|- The samples are 
identified in Table The Bai_a;Ca2:F2 samples showed 
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FIG. 1: Raman and x-ray determination of the mole per- 
cent CaF2 composition in the mixed crystals. Open points: 
Bai_a;Caa;F2; filled points: Pbi_iCaiF2. Circles (squares) 
are from Raman T25 measurements along the sides (ends) of 
the samples. Diamond points are from powder x-ray scatter- 
ing of end-cut pieces. 



TABLE I: Nominal and actual dopant concentrations XNom. 
and X, in mole % for individual mixed crystals; and the di- 
electric constant eo for each host. 



eo = 6.94 


^Nom. 


X 


eo = 26.3 


^Nom. 


X 


BaCal 


2 


2.2 ±0.5 


PbCal 


2 


2.5 ±0.3 


BaCa2 


4 


3.1 ±1.1 


PbCa2 


10 


3.9 ±0.9 


BaCaS 


6 


5.4±1.1 


PbCa3 


6 


4.0 ±1.0 


BaCa4 


8 


5.4 ±1.2 


PbCa4 


4 


4.2 ±0.2 








PbCa5 


8 


5.9±1.6 



significant variation over the sample length for greater 
than 2 % CaF2, consistent with the onset of the misci- 
bility gap [10]; while the Pbi_xCaxF2 samples were in- 
homogeneous above about 4 % CaF2, again consistent 
with the phase dia gram showing segregation starting at 
4-6 mol % CaF2 [llj. The mean CaF2 concentrations 
were extracted by averaging the Raman data taken along 
the sides of each sample (Table IJ). Laue back-scattering 
measurements indicated that the samples are composed 
of few-mm single crystal domains with slight reorienta- 
tions between neighboring crystals. 

Typical low temperature (1.5 K) THz absorption spec- 
tra over a decade in frequency are plotted in Fig. |21 for 
some of these mixed crystal samples. These spectra show 
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FIG. 2: THz absorption coefficient of mixed crystals versus 
frequency at 4.2 K. No sharp defect-induced absorption peaks 
are observed. Sample concentrations are given in Table 



that no sharp impurity-induced absorption lines appear 
from 2-25 cni^^, including the Bao.996Cao.oo4F2 sample. 
Such sharp features are typical in the low concentration 
spectra of aliovalent mixed crystals 

When TLS make a significant contribution to the ab- 
sorption, then increasing the sample temperature pro- 
duces the characteristic bleaching observed in glasses 14] . 
The predicted temperature-induced change in the ab- 
sorption coefficient a, is 



Aa(Tji,r) 



4^^ 



wP(^) fj^lito) Ap{Tr,T), 



V 3 

(1) 

where Ap{Tr,T) = [ta,iih{huj/2kTR) ~ ta,iih{hLu/2kT)]. 
The quantity P{uj)nl{u!) is the optical TLS density-of- 
states (ODOS). Its frequency dependence is due to a 
cutoff in the high frequency tunneling energy and the 
onset of excited state transitions |^. Taking the limit 
P{uj 0) yields P{0), the constant spectral density of 
states per unit frequency range, is the electric dipole 
transition moment, c is the speed of light, the term in 
parentheses is the local field correction, eg is the dielec- 
tric constant, and k is Boltzmann's constant. 

The temperature dependent spectra between 1.5 and 
10 K are used in conjunction with Eq. ^to extract the 
ODOS P{uj)^'^{uj), which is graphed versus frequency for 
both series in Fig. |3 The Bai_a;CasF2 spectra show 
two-level behavior increasing monotonically with concen- 
tration above a certain minimum value. For the lower 
concentration BaCal and BaCa2 samples, the TLS dis- 
tribution is not yet completely glass-like since it does 
not extend below about 4 cm~^. Nonetheless, there is 
a significant range of two-level behavior and these data 
are used in Fig. 0] For the BaCa3 and BaCa4 samples, 
broad distributions of TLS are already present, although 
saturation has not been reached. The Pbi_a;Ca2:F2 data 
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FIG. 3: Measured TLS optical density of states of the mixed 
crystals versus frequency. The data are obtained from sin- 
gle parameter fits to the temperature dependence of the THz 
absorption coefficient, (a) 'Q&i-xCa.x^2 and (b) 'Phi-xCa.x^i- 



shown in Fig. Ol^b) demonstrate negligible ODOS for 
PbCal (2.5 %) and a non-monotonic concentration de- 
pendence. For the Pbi_a;Caa;F2 samples, the onset of 
TLS behavior is very sharp and therefore sensitive to the 
precise sample concentrations. The maximum measured 
ODOS for Bai_2;Caa;F2 is about one tenth of that in 
saturated mixed crystals and glasses 0, and the ODOS 
in Pbi_a;Caa;F2 is about one hundredth of the saturated 
values. 

Studies with doped alkali halide crystals have already 
demonstrated that for a given lattice containing a smaller 
size dopant cation that although single ions may have in- 
sufficient volume to tunnel, pairs of such ions can [isj . 
In both of the systems studied here the dopant ion is the 
smaller; but there is no evidence of single ion or pair tun- 
neling. We suggest that the TLS effect observed here is 
associated with a local statistical concentration of dopant 
ions exceeding a critical value Xc- 

Since the probability that a given dopant ion is lo- 
cated in a specific subvolume is (subvolume/ crystal vol- 
ume), the characteristic subvolume for the disordered 
state needs to be identified. The first step is to intro- 
duce a Voronoi tessellation centered on all dopant sites 
that fills all space |0| . In this way two relevant volumes 
are defined; Vi is the Voronoi polyhedron volume at site 
I, and V is the expectation value of Vi. The normalized 
Voronoi volume at site i is then Zi = Vi/V. Assuming 
random dopant locations in the crystal and neglecting 



TABLE IL Critical concentration model fit parameters. 



Crystals 


[PM?].at(xlO-5) 


Xc (mol %) 


(CaF2)i-,(LaF3), 


113±4 


0.56 ±0.06 


(SrF2)i-4LaF3):. 


56 ±6 


10 ±2 


(BaF2)i_a;(LaF3).: 


570 ± 60 


9.5 ±1.4 




1.6 ±0.3 


3.2 ±0.3 


Pbi_,Ca,F2 


1.6 


7.5 ±0.05 



lattice discreteness, a Poisson Voronoi diagram follows. 
One feature of the Poisson diagram is that the proba- 
bility density distribution for Zi can be approximated by 
the gamma distribution 



P{Zi) 



b^T (a) 



(2) 



where in 3-D a — 5.6333 and h — 0.1782. Summing all 
Voronoi subvolumes less than a fixed critical value 14, 
which is the fitting parameter, and defining Zc — Vc/V 
gives the probability W of locally exceeding a critical 
concentration 



J dz,z,p{z,). (3) 



In order to make contact with the ODOS, the satura- 
tion level [Pii^]sat is introduced as a second parameter 



so that the final expression is 



[P^^b\so.f I dZiZip{Zi). 



(4) 



The critical concentration can be computed from the fit 
parameter Vc by Xc = VhostlVc, where Vhost is the volume 
of the host crystal primitive cell. 

This critical concentration model is quite general and 
may be tested against all the THz data on TLS in mixed 
crystals. Fits have been carried out for three aliovalent 
impurity systems, namely CaF2, BaF2, and SrF2 doped 
with LaFs and for the two isovalent systems studied here. 
The results are shown in Fig. ^and the parameters given 
in Table ITU The model reproduced well both the delayed 
increase in TLS absorption with concentration as well 
as the saturation at higher concentration. The appar- 
ent underestimation by the model at low concentrations 
for CaF2:LaF3 may arise because narrow defect-induced 
tunneling states appear simultaneously with the incipi- 
ent broad TLS absorption. In addition, clustering effects 
would result in higher than random appearance of regions 
with concentrations exceeding the critical value. 

For the isovalent impurities studied here, the critical 
concentration model also shows qualitative agreement 
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FIG. 4: Critical concentration model fits to the TLS optical 
density of states for both aliovaleut and isovalent mixed crys- 
tals. Data points are a weighted average from 3-10 cm~^. 
The samples are identified in the figure. 



with the experimental data. The saturation level ob- 
tained by fitting the Bai_a;Caj:F2 data is two orders of 
magnitude less than for the aliovalcnt mixed crystals. 
Most probably this is the result of a smaller dipole mo- 
ment Hb since the absorption by TLS created in isovalent 
systems is not enhanced by charge separation. The large 
uncertainty in the concentration of Pbi_a;Caa;F2 samples 
did not permit a reliable determination of [P^iDsat so the 
Bai_a;Ca2:F2 value was applied. 

It has been shown that two isovalent mixed crystals, 
Bai_a;Caa;F2 and Pbi_a;Caa;F2 do contain significant al- 
though not saturated TLS optical density of states. The 
delayed but rapid onset of TLS behavior with increas- 
ing dopant concentration in these minimally disordered 
monatomic substitutional systems has been successfully 
modeled with a critical concentration parameter related 
to the statistics of the dopant distribution in the lattice. 
The other parameter, the saturation value, indicates that 
only a small fraction (~ 1/1000) of the critical concen- 
tration regions contribute to the collective tunneling phe- 
nomenon. This fact suggests a complex energy landscape 
even for these simple lattice systems. 



Finally, the observation of a step-like appearance of 
TLS with the number density, independent of the local 
concentration as soon as it exceeds the threshold value, 
implies that there are no intermediate states either TLS 
are present at the saturation level or they are absent. 
The success of this very simple model for mixed crystals 
points to a likely connection with the universality of TLS 
found in glasses. 
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